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In this paper we introduce Parity-Time (VT) symmetric perturbation to a one-dimensional Lieb 
lattice, which is otherwise P-symmetric and has a flat band. In the flat band there are a multitude 
of degenerate dark states, and the degeneracy N increases with the system size. We show that the 
degeneracy in the flat band is completely lifted due to the non-Hermitian perturbation in general, 
but it is partially maintained with the half-gain-half-loss perturbation and its “V” variant that we 
consider. With these perturbations, we show that both randomly positioned states and pinned 
states at the symmetry plane in the flat band can undergo thresholdless VT breaking. They are 
distinguished by their different rates of acquiring non-Hermicity as the PT-symmetric perturbation 
grows, which are insensitive to the system size. Using a degenerate perturbation theory, we derive 
analytically the rate for the pinned states, whose spatial profiles are also insensitive to the system 
size. Finally, we find that the presence of weak disorder has a strong effect on modes in the dispersive 
bands but not on those in the flat band. The latter respond in completely different ways to the 
growing PT-symmetric perturbation, depending on whether they are randomly positioned or pinned. 


I. INTRODUCTION 

Systems that exhibit flat bands have attracted consid¬ 
erable interest in the past few years, including optical 
mm and photonic lattices lam, graphene 13 0 , su¬ 
perconductors mm, fractional quantum Hall systems 
prong and exciton-polariton condensates min]- One 
interesting consequence of a flat band is the different 
scaling properties of its localization length when com¬ 
pared with a dispersive band jT8H20] . due to a multitude 
of degenerate states in the flat band. This degeneracy 
has another important implication on Parity-Time (VT) 
symmetry breaking [2TH39] : it was found recently that 
the degeneracy of the underlying Hermitian spectrum, 
before any T-breaking perturbations are introduced, de¬ 
termines whether thresholdless VT symmetry breaking 
is possible [40] . Therefore, it is interesting to probe the 
interplay of the large degeneracy in a flat band system 
and VT symmetry breaking. In particular, we would 
like to know whether the degeneracy is completely lifted 
due to a PT-symmetric perturbation, and whether VT 
symmetry breaking depends on the evenness (oddness) of 
the degeneracy and its magnitude that grows with the 
system size. In addition, because a flat band makes the 
underlying Hermitian system more susceptible to disorder, 
it is equally important to understand the role of disorder 
on VT symmetry breaking in a flat band system. 

Using a quasi-one-dimensional (quasi-ID) Lieb lattice 
(see Fig. [l]) and a half-gain-half-loss perturbation (includ¬ 
ing its variant, the “V” configuration to be introduced 
below), we show in this paper that two different sce¬ 
narios of thresholdless VT symmetry breaking can take 
place, depending on whether the degeneracy N in the 
flat band is even or odd. When N is odd, all but one 
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flat band modes enter the PT-broken phase at the in¬ 
finitesimal strength of the PT-symmetric perturbation. 
They form two (N — l)/2 degenerate branches, each con¬ 
fined strictly to either the gain half or the loss half of 
the lattice, and their non-Hermicity equals the strength 
of the PT-symmetric perturbation, which we denote by 
r. The remaining flat band mode experiences an excep¬ 
tional point of order 3 [41] . via the coupling with two 
dispersive band modes. When N is even, all flat band 
modes experience thresholdless VT breaking, but now 
they exhibit four branches. Two branches are (TV/2 — 1) 
degenerate and have the same properties as those in the 
TV-odd case. There is only one mode in each of the re¬ 
maining two branches, and they form a PT-symmetric 
doublet, pinned at the symmetry plane (x = 0) with 
exponential tails in both the gain and the loss halves. 
Surprisingly, we find that this localization is not due to 
the half-gain-half-loss nature of the PT-symmetric pertur¬ 
bation as previously found [42j , but rather a result of the 
point defect at x = 0 of this perturbation to satisfy the 
VT symmetry. Therefore, this spatial profile is insensitive 
to both r and the system size when N is large, and it 
can be reproduced in the underlying Hermitian system at 
r = 0. Using a degenerate perturbation theory, we derive 
this localization length and the rate these doublet states 
acquire non-Hermicity. 

Finally, we show that the presence of weak disorder has 
a strong effect on modes in the dispersive bands but not on 
those in the flat band. The finite PT-transition thresholds 
of the former are smoothed out, while the thresholdless 
VT breaking of the latter is largely preserved. In addition, 
the flat band modes respond in completely different ways 
to increasing PT-symmetric perturbation when there is 
disorder. The PT-symmetric perturbation has little effect 
on them if they are already localized to either half of the 
lattice. Otherwise the perturbation forces them to pick 
a side, unless they are the doublet states, which evolve 
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FIG. 1. (Color online) Inset: A quasi-ID Lieb lattice with 3 
overlapping “ A-B-C-B-A ” sublattices. Two dark states are 
shown schematically by circular clouds, which also serve as an 
example of a RT-symmetric “V” perturbation to be discussed 
in Section II. Main figure: Energy eigenvalues in a Lieb lattice 
with 10-fold degeneracy in the flat band. There is one more 
mode in each of the two dispersive bands. In all figures the 
energy is scaled by the coupling G(= J), where the onsite 
energies are all shifted to zero (uja — ujb — wc = 0). 


from Anderson localized states [43] to those pinned by 
the point defect at x = 0. We note that different from 
Refs. 2F 4T 45], the system we consider here has a flat 
band before the VT -symmetric perturbation is introduced, 
instead of being the result of PT-symmetry breaking. 

II. THRESHOLDLESS VT BREAKING IN A 
LIEB LATTICE 

In two-dimensional systems of a finite size, the maxi¬ 
mum degeneracy due to a point group is 2. In a flat band, 
however, the degeneracy can be arbitrarily large, because 
it increases with the system size. Take a quasi-ID Lieb 
lattice for example (see Fig. [lj, where an upright A site 
is decorated to every other site of a ID chain (i.e., the 
B sites). A flat band exists when the onsite energy of A 
and C sites are the same, uja = uc = {h = 1), and 
cjo gives the flat band energy. It is gapped from the two 
dispersive bands by G, the real-valued coupling between 
two neighboring A and B sites. For simplicity we take 
ujb = cjo as we ll- Because cjo shifts all eigenvalues by the 
same amount and has no effect on the eigenstates, we 
take it to be zero. The degeneracy N of the flat band in 
this Lieb lattice equals the number of existing dark states, 
which we will denote by V(x) due to their geometry (see 
the inset in Fig. [l]). Each dark state has a nonvanishing 
amplitude only at a C site and two nearest A sites: 


where the five elements are the amplitudes of the wave 
function on a “ A-B-C-B-A ” sublattice, to which we 
will refer as an V sublattice below. The B sites are 
black in V(x) because the tunneling probabilities from 
the neighboring A and C sites cancel each other. The 
argument x in V(x) is the position of the central C site 
of an U sublattice in the units of the lattice constant. It 
takes integer or half integer values, depending on whether 
N is odd or even, and the center of the lattice is placed 
at x = 0. We note that the number of (overlapping) U 
sublattices equals the number of dark states, and each 
end of the lattice is terminated on A and B sites. The 
superscript “T” in Eq. © denotes the matrix transpose, 
and J is the real-valued coupling between two neighboring 
B and C sites. Below we drop the vector symbol of V(x) 
without causing ambiguity. 

The tight-binding Hamiltonian of the Lieb lattice can 
be written as 

Ho =J2^iZj)(Zj\ + £ [G\Aj){Bj\ + h.c.} + 

j 

E [ J ( \ B i){Cj\ + \Cj)(B j+ i\) + h.c.], (2) 

j 

where j runs through all unit cells, Z m A, B, G, and h.c. 
denotes Hermitian conjugate of the other terms in the 
square brackets. We then introduce the non-Hermit ian 
perturbation irHi, which is PT-symmetric about x = 0. 
r is the overall strength of the perturbation, and Hi is 
a diagonal matrix with positive (gain), negative (loss), 
and zero (no gain or loss) elements. To satisfy the VT- 
symmetry, no gain or loss is introduced on the lattice 
sites right at the symmetry plane at x » 0, which include 
one A and one B site when N is even and a single G site 
when N is odd. 

The degeneracy of the flat band modes is completely 
lifted by Hi in general as they undergo thresholdless 
VT symmetry breaking, independent of whether N is 
even or odd (see Fig. [2]). The simplest PT-symmetric 
perturbation, with uniform gain on one half of the whole 
lattice and the same amount of loss on the other half, in 
fact partially maintains the degeneracy of the flat band. 
A variant of this half-gain-half-loss configuration with the 



FIG. 2. Thresholdless VT breaking in a Lieb lattice that 
completely lifts the degeneracy in the flat band with a general 
RT-symmetric configuration for (a) iV = 5 and (b) IV = 6. 


V(x) = [J, 0 , -G, 0 , J] T 


( 1 ) 
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FIG. 3. (Color online) Thresholdless VT breaking in a Lieb 
lattice with the “V” configuration for (a) TV — 5 and (b) 
TV — 6 . (a) 4 flat band modes enter the PT-broken phase at 
r as 0 (thick solid lines). The other flat band mode and all 12 
dispersive band modes remain in the symmetric phase in the 
range of r shown (thin solid line), (b) All 6 flat band modes 
enter the PT-broken phase at r = 0, 4 with k± = dbr (thick 
solid lines) and 2 with 0 < \k±\ < r (dashed lines). All 14 
dispersive band modes remain in the symmetric phase in the 
range of r shown (thin solid line). 


same property is the “V” configuration (see Fig. [lj inset), 
which differs by having no gain or loss on all the B sites 
(where the flat band modes are dark). 

To see how the half-gain-half-loss configuration and the 
“V” configuration partially maintain the degeneracy of the 
flat band, we first note that a pair of V(x), V(—x) modes 
are VT- symmetric partners, i.e., V(x) = VTV(-x), be¬ 
cause they are real-valued wave functions (see Eq. 0). In 
addition, they are also eigenfunctions of Hi for these two 
configurations, with opposite eigenvalues (± 1 ) when they 
do not overlap with the A and C site at x = 0. There¬ 
fore, they undergo thresholdless VT symmetry breaking 
simultaneously as r increases from 0 , and the imaginary 
parts of the corresponding energy eigenvalues are nothing 
but k± m Im[d)±] = ±t (see Fig. [3]), independent of the 
system size and the couplings G, J. In other words, their 
non-Hermiticy (given by \k±\) is the same as the strength 
of the PT-symmetric perturbation (i.e., r). There are 
(.TV — l)/2 such pairs of flat band modes when TV is odd 
(e.g., 2 and 3 pairs for the N m 5,7 cases shown in 
Figs. [3^a) andj^a)) and (TV/2 — 1) such pairs when TV 
is even (e.g., 2 and 3 pairs for the TV = 6 , 8 cases shown 
in Fig. [3^b) and [ 6 ^b)), which all behave in the same 
way. As a result, each of the two k± = ±r branches has 
(TV — l )/2 (TV-odd) or N/2 — 1 (TV even) degeneracy in 
the VT- broken phase. Such a large number of degenerate 
states in the PT-broken phase has not been reported 
before. 

For the single dark state V(x «■= 0) that is left alone 
in the TV-odd case, it is no longer an eigenstate of Hi 
since there is no perturbation on the G site at x = 0 
in order to respect the VT symmetry. Hence it needs 
to couple to modes in the dispersive bands to break the 
VT symmetry, which then leads to a finite threshold 
in terms of r. In fact at this PT-symmetry breaking 
threshold lies an exceptional point of order 3 (EP 3 ) |41j , 
where the energy eigenvalues and wave functions of three 
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FIG. 4. (Color online) EP 3 in a Lieb lattice with the “V” 
configuration and N — 7. The real and imaginary parts of 
the energy eigenvalues of two dispersive band modes and the 
solitary dark state are shown in (a) and (b) respectively, with 
the EP 3 marked by the black dot. False color plots of the 
absolute value of their wave functions at (c) r = 0 and (d) 1.58 
are shown. The white U sublattice in (c) indicates different 
lattice sites. 


modes coalesce, including two dispersive band modes 
and the solitary dark state (see Fig. [4]). Nevertheless, 
a hybridized mode remains in the PT-symmetric phase 
after the EP 3 , as if the solitary dark state remained in 
the PP-symmetric phase. We note that this EP 3 does 
not occur if the onsite energy ujb is detuned from uja and 
l dc or with the half-gain-half-loss configuration. 

The two dark states V(x = — 1 / 2 ) and V(x = 1/2) in 
the TV-even case are not eigenstates of Hi either, due to 
the absence of perturbation on the A site at x = 0 in the 
“V” configuration (and the B site at x = 0 as well in the 
half-gain-half-loss configuration) in order to respect the 
VT symmetry. We will refer to these two states as the 
doublet states. Although they do not belong to the (TV /2 — 
1 ) degenerate manifolds of the two k± = ±r branches, 
they still experience a thresholdless VT breaking, but 
their \k,±\ is smaller than r (see Fig. |3jb)). Interestingly, 
we find that the slope of these \k±\ is also insensitive to 
the system size TV, when TV > 2. More specifically, this 
slope is about 0.708, 0.675, 0.671, 0.670 when TV = 2,4, 6 , 8 , 


measured by k'+{t = 0 . 1 ) ~ ^re¬ 
configuration (see Fig. |3^b)), which 
we will focus on below. 


with the “V” 

T = 0.1 

is the configuration 


To understand this property, we resort to two ap¬ 
proaches, a qualitative one based on directly visualizing 
the spatial profiles of the doublets states, and a quantita¬ 
tive one based on a degenerate perturbation theory. As 
Fig. J5^a) shows, the doublet states are pinned near x = 0 
and have tails in both the gain and the loss halves of the 
lattice. We find that these spatial profiles are also insensi¬ 
tive to the system size when TV is large enough (> 2; see 
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FIG. 5. (Color online) Doublet states in the PT-broken phase 
with the “V” configuration, (a) False color plot of the absolute 
value of their wave functions for N = 6. The vertical dashed 
line marks x = 0, and the color scale is the same as in Figs.[4|c) 
and (d). (b) Same as (a) but for N = 20. (c) Exponential tails 
of the amplifying doublet state shown in (b) (solid lines). Only 
the absolute values of the wave function on the C sites are 
shown. Circular dots and open triangles are calculated at r — 
0.002 and 0.2, respectively. The solid lines are proportional to 
exp(— |x|/1.03) shown for comparison. The dashed line shows 
the defect state at r = 0 by increasing uja{x — 0) to 0.1. 


Figs.^a) and (b)), which explains the same property of 
their k±. In fact, we find the tails of their spatial profiles 
decay exponentially away from x = 0 (see Fig. [5jc)). One 
might expect that these exponential tails are a result of 
gain and loss, as it has been shown that the wave function 
in the PT-broken phase tends to peak at the gain and 
loss interface [42] . However, we find that the exponen¬ 
tial tails are largely independent of the PT-symmetric 
perturbation strength, either. As Fig. [5^c) shows, the 
tails are captured well by an exponent of — |x|/1.03, and 
there is little change in the wave function when r changes 
from 0.002 to 0 . 2 . Instead, we find that this localization 
length of £ = 1.03 is due to the point defect of the “V” 
configuration at x = 0 , which has neither gain or loss on 
the A site at x = 0 in order to respect the PT symmetry 
as mentioned. If we introduce a point defect at x = 0 in 
the Hermitian system at r = 0, for example, by increasing 
uj a here to 0 . 1 , we recover the same localization length 
(see the dashed line in Fig. [5jc)). 

Before we derive the perturbation theory for the doublet 
states, we briefly discuss the modes in the dispersive bands. 
They have high PT-transition thresholds when N is small. 
This observation can be understood in the following way: 
their eigenvalues are distributed in \u\ G [cjo + G,cj o + 
VG 2 + 4 J 2 ] (see Fig. [lj and are well separated from each 
other when N is small; their thresholds are proportional to 
these spacings, in the simplest case of two-mode coupling 
m- As N increases, their energy spacings reduce and 
so do their PT-transition thresholds. For example, no 
dispersive band modes enter the PT-broken phase in 
Fig-! but their lowest threshold reduces to r ^ 0.090 
and 0.072 for N = 7 and 8 , respectively (see Fig. [ 6 ]). Due 
to the symmetry of the dispersive bands about cjo = 0 
(see Fig. [l]), when one pair of dispersive band modes enter 



FIG. 6. (Color online) Same as Fig. [3] but with (a) N = 7 
and (b) N = 8. (a) 4 and (b) 8 dispersive band modes enter 
the PT-broken phase before r — 0.1 (dash-dotted lines). Note 
that the features of the flat band modes (solid and dashed 
lines) are the same as in Fig. [3] The horizontal thin lines 
represent modes in the PT-symmetric phase. 


the PT-broken phase at some finite value of r, there is 
always another pair that do the same at exactly the same 
r but on the opposite side of the flat band energy ujq; 
their k± are nevertheless the same. Different from the 
EP 3 scenario depicted in Fig. [4] in which two dispersive 
band modes are also involved, here the real parts of their 
energy eigenvalues are different from ujq in the PT-broken 
phase (not shown). 


III. DEGENERATE PERTURBATION THEORY 

Now we turn to the perturbation theory for a more 
quantitative understanding of the doublet states, and 
hence the N considered in this section is even. We start 
by diagonalizing the perturbation Hi in the TV-fold flat 
band subspace of Ho, which is required in the degenerate 
perturbation theory [46] . In the N = 2 case, we define 
the basis in this subspace as (j>^ = V(—1 /2) + aV( 1/2) 
and 02 = bV(—1 /2) + k(l/2), where we can take a,b 
to be real numbers, because H 1 (without the factor ir) 
is Hermitian. The superscript enclosed in parentheses 

indicates the value of TV. The above requirement then 
(2) (2) 

means ( 0 ^ ; |iii |02 / = 0 , and an additional constraint is 
(0i 002 )= 0. We note that this requirement is different 

from demanding that 0 ^ , 02 are the eigenstates of H 1 ; 
the latter is sufficient but unnecessary. 

Since the flat band states are dark on B sites, we can 
project the wave functions and Hi to the Hilbert space 
of A and G sites, leading to V(— 1 / 2 ) = [J, —G, J, 0 , 0 ] T , 
V"(1/2) = [0,0, J, —G, J] T , and Hi = diag(l, 1,0, —1, —1) 
when TV = 2, where “diag” stands for a diagonal matrix. 
Using (0^ 2 ^ |-H"i 1 02 ^) = 0 we find a = b, and (0^ 2 ^ 1 02 2 ^) = 
0 leads to 



Note that the two values of a are reciprocal to each other, 

( 2 ) 

which means that if one leads to (j >\ ) , the other one gives 
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( 2 ) 

02 . By neglecting the coupling to the dispersive band 
modes, we then derive 


.( 2 ) 






<d 2 ) id 2) > 


a 2 - 1 


2 a 


J 2 


3 J 2 + G 2 


(4) 

for the doublet states, which gives a slope of 0.707 when 
G = J = 1. Using 0^ i n the expression above leads to 
the same result. Previously we have mentioned that the 
numerically obtained slope in this case is 0.708, calcu¬ 
lated using k/,(t = 0 . 1 ) ~ 


. If we reduce the 


T = 0.1 

perturbation strength at which the slope is calculated, for 

example, to r = 0 . 02 , we find = 0.707, which 

T T = 0 - 02 _ 

agrees nicely with the analytical result given by Eq. Q . 

Since the slope of \k±\ given by Eq. Q is independent of 
r, we know that any change to the slope must be a result 
of coupling to the dispersive bands, which is neglected 
in the derivation of Eq. Nevertheless, the small dif¬ 
ference between a t r = 0.02 and 0.1 indicates that 

T 

such coupling is weak. 

For TV > 2, the construction of the first TV — 2 basis 
functions in the TV-fold flat band subspace is straightfor¬ 
ward: any dark states that do not overlap with x = 0 are 
eigenfunctions of Hi, as we have mentioned. These TV — 2 
dark states just need to be linearly superposed properly 
such that they form TV — 2 orthogonal basis functions, 
and there is more than one way to achieve it. For the 
remaining two doublet states, we use mathematical in¬ 
duction to find their approximate forms. Assuming that 
we have solved the (TV — 2) case and found the correct 
doublet states <p[ N ^, 0 ^ ^ ? we then approximate the 
correct doublet states in the TV case by 


« J>[ N -V + cV{-N/ 2) + dV{N/2), (5) 

4 N) * 4 N ~ 2) + dV(-N/ 2) + cV(N/ 2), ( 6 ) 

where c, d are two real numbers. The basic assump¬ 
tion is that the central part of (j)[ N ^ , 0 ^ near x = 0 
is insensitive to TV, as we have seen in Figs. [5ja) and 
(b). It is easy to check that (0^ |iii 10^) = 0, using 

{(f)[ N ~ 2 ' , \H 1 \4 N ~ 2 ' > ) = o, and by requiring (4V IfiV) = 
0 , we find 

cd(2J 2 + G 2 ) + 2 J(c 5 {n ~ 2) + de {N ~ 2) ) = 0, (7) 

where e^ N ~ 2 \S^ N ~ 2 ^ are the first (last) and last (first) 
elements of <p[ N ~ 2 ^ ( 02 /V_2 ^). An additional constraint 
imposed by the degenerate perturbation theory is that 
(f )[ N ^, 0^ need to be orthogonal to the first TV — 2 ba¬ 
sis functions in the TV-fold flat band subspace. Let us 
pick, for example, one such basis function simply as 
03 ^ = U(—TV/ 2 ) (which is an eigenfunction of -Hi), and 
we immediately find d = —T N ~ 2 ^/3 from ( 0 ^ | 0 ^) = 0 

and c — —e^ N ~ 2 ^/3 from ( 03 ^ 102 ^) = 0 when G 
J 1 . We note that we arrive at the values of c 


and d with other choices of 0^, such as V(N/2) or 
V(— TV/2) — V(N/2). Once c and d are found, we can 
immediately construct the leftmost and rightmost ele¬ 
ments of and 0 ^ to derive the localization length 
for the doublet states. For example, the first four ele¬ 
ments of 0^ are cJ, —cG,cJ-\-e^ N ~ 2 \ —e^ N ~ 2 ^G/J. The 
localization length £ can then be defined as 


r 1 = In 


e (JV- 2 ) 

cJ 


(8) 


using the values of the wave function on the first C site 
(i.e., — cG ) and the second C site (i.e., —e^ N ~ 2 ^G/J) from 
the left. Equation ^ gives £ = In 3 ~ 1.10, which agrees 
reasonably well with the result of numerical fitting of 
the exponential tails shown in Fig. |5^c), i.e., £ = 1.03. 
In the Appendix we give another way to estimate this 
localization length. 

For the slope k±(t) of the doublet states, a recursive 

relation can then be formulated. Assuming 0^ is ampli¬ 
fying, we find 

+ ^rvr 1 ) 

O ,(JV-2) 5 (iV-2)2_ e (]V-2)2 

_ + (4 w - 2 y- 2) > rq x 

£(iV-2)2 +e (jV-2)2 • V*) 

<^- 2 ) l^- 2) > 


It is clear that (and similarly becomes insen¬ 

sitive to TV as TV becomes large, because e^ N ~ 2 \ 5^ N ~ 2 ^ 
are at the very ends of the exponentially decaying tails 
of 0^- 2) , i.e., |y(iV-2)| 2 ,| e (iV-2)|2 <<c < 0 ( iV - 2 )| 0 (^- 2 ) > , 

leading to a converging series w Using 

= 0.707 given by Eq. (14), we successively find 

= 0.673, = 0.669, = 0.669, which agree 

well with the previously mentioned numerical values. 


IV. EFFECT OF DISORDER 

Unlike the doublet states, the VT -broken dark states 
in the two degenerate manifolds can be randomly posi¬ 
tioned along the lattice, but they are strictly confined 
on either the gain or the loss half of the lattice. Due to 
their degeneracy, any arbitrary superpositions of them 
are still eigenstates of the system, and one can always 
find pairs of them such that they are VT -symmetric part¬ 
ners of each other. However, the relevant superpositions 
are determined by the disorder in the system. Here we 
consider weak disorders by including a white noise with a 
uniform distribution in [— W/2, W/2] (W <C G, J) around 
cjo o n each lattice site. As the lower mode in Figs. [7] 
(a)-(c) shows, if a flat band mode in the two degenerate 
manifolds is already confined to the loss (or gain) half of 
the lattice at small r, increasing r has little effect on its 
spatial profile. However, if at small r one of these modes 
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FIG. 7. (Color online) (a)-(c) Evolution of the dark states in 
the N/ 2—1 degenerate flat band manifolds from r — 5 x 10 -4 , 
0.01, to 0.075. Two typical examples are shown, (d)-(f) Same 
as in (a)-(c) for the doublet states. The color scale and legend 
are the same as in Figs.[5ja) and (b). The amplitude of onsite 
white noise is W — 0.05 and N — 20. 


occurs in both the gain and loss regions, it will be forced 
to take either the gain half or the loss half, as the upper 
mode in Figs. [7] (a)-(c) shows. Here the prevailing effect 
of irH\ is its non-Hermicity, separating the lattice into a 
gain half and loss half. 

In contrast, the prevailing effect of irH\ on the dou¬ 
blet states is its point defect mentioned previously. At 
very small r, the doublet states are localized by disorder, 
i.e., they are Anderson localized m, and the localization 
length can be rather long in a particular disorder realiza¬ 
tion (although on average the localization length is not 
much longer than that determined by the point defect 
ED]). As r increases, the point defect in Hi gradually 
prevails over the white noise disorder, and the doublet 
states approach their spatial patterns in the absence of 
FF, as can be seen from Figs. 0e) and[7p) in comparison 
with Fig. [5] In other words, the VT- symmetric relation 
0 i = VT (/)2 of the doublet states is restored as r increases. 
The same cannot be said about the other flat band modes 
in general, due to the presence of disorder that breaks 
the VT symmetry. 



FIG. 8. (Color online) Same as Fig. [6] but with disorder. The 
amplitude of onsite white noise is W — 0.05. 


Finally, we find that weak disorder has a stronger effect 
on the dispersive bands than the flat band in terms of 
VT breaking. As Fig. [ 8 ] shows, the disorder smoothes out 
the PT-transitions of the dispersive band modes in Fig. [ 6 ] 
but the thresholdless VT breaking of the flat band modes 
is largely unaffected. These contrasting behaviors are 
due to the different routes to PT-symmetry breaking, i.e., 
whether an exceptional point is involved. On the one hand, 
an exceptional point is very sensitive to perturbations, 
and it quickly becomes an avoided crossing in the complex 
eigenvalue plane when the PT-symmetry of the system 
is lifted by the weak disorder. In fact, each bifurcation 
of Im[cj] we have seen in Fig. [ 6 ] from a finite r is for 
two pairs of dispersive band modes as mentioned at the 
end of Section [II] Now in the presence of weak disorder, 
the sensitivity of the exceptional points breaks each of 
these bifurcations into two, and we see four and eight 
trajectories of Im[ cj] in Figs. [ 8 ^a) and (b) (dash-dotted 
lines) instead of two and four in Figs.^a) and (b). On 
the other hand, the thresholdless VT -symmetry breaking 
of the flat band modes does not involve an exceptional 
point; its role is replaced by the Hermitian degeneracy 
at r = 0 [40] . In this case there is no singularity, and 
we have seen that the first-order perturbation theory 
describes the thresholdless PT-breaking well in Section 


III Now with the weak disorder, its leading (i.e., linear) 


effect is merely perturbing the real part of the energy 
eigenvalues, similar to how the PT-perturbation varies 
the imaginary part of the energy eigenvalues. Therefore, 
these two perturbations are independent processes to 
the leading order, and the thresholdless PT-breaking is 
largely unchanged with the weak disorder. As the strength 
of PT-breaking increases with r, the effect of the weak 
disorder is further suppressed. This can be seen from 
the following observation. We find that one flat band 
mode with ~ —r in Fig. [ 8 jb) deviates slightly from 
the other 2 in this branch at small r but rejoins the latter 
as r increases. This behavior is due to its spatial profile 
change caused by the increasing VT perturbation, similar 
to that of the upper mode in Figs. [7](a)-(c): when this 
mode has a small portion in the gain half of the lattice 
at small r, overall it experiences less loss and hence its 
negative is larger than —r. As its spatial profile is 
squeezed into the loss half of the lattice by the growing 
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VT- symmetric perturbation, it experiences more loss than 
before and its /c_ approaches its minimum value —r. We 
can also infer that the other 2 modes in this branch are 
localized already in the loss half when r is small, which is 
verified by inspecting their spatial profiles (not shown). 


V. CONCLUSION 

We have mentioned that the degeneracy of the flat 
band modes is completely lifted by a VT- symmetric per¬ 
turbation in general as they undergo thresholdless VT 
symmetry breaking. An exception takes place in the half- 
gain-half-loss configuration and its “V” variant, where 
two different scenarios of thresholdless VT symmetry 
breaking occur in the flat band, depending on whether 
the degeneracy N in the flat band is even or odd. The 
two degenerate manifolds with k± = =br always exist and 
undergo thresholdless VT breaking. While their center 
positions are random, they are confined to either the 
gain half or the loss half of the lattice in the absence of 
disorder. This feature holds even with weak disorders, 
when the PT-symmetric perturbation is strong enough. 
In contrast, the PT-symmetric doublet states only exist 
when N is even, and they display weaker non-Hermicity 
than the degenerate manifolds. These doublet states are 
pinned at the symmetry plane with exponential tails in 
the gain and the loss halves. These tails are the result 
of a point defect in the PT-symmetric perturbation at 
x = 0, instead of its half-gain-half-loss nature as previ¬ 
ously found [42]. Weak disorders may disturb their spatial 
profiles at small PT-symmetric perturbation, but as the 
latter increases, this feature is restored, together with the 
PT-symmetry relation between the two doublet states. 


where H' represents any perturbation to the diagonal 
elements of Hq (including the “V” configuration and the 
Hermitian defect discussed in Fig. [5]) and Zj(Z = A , P, C ) 
is the wave function on the j th site of type Z. A zj = 
uj — cozj , where uj is an eigenvalue of H and ujzj is the 
onsite energy of Zj. ujzj can have an imaginary part to 
represent gain or loss if present. 


By solving the Eqs. (10) and ( [IT] ), we find 


B 0 = 


ACj + Cj- 1 ) 

A G 2 ' 


(13) 


By substituting Bj and Pj+i in Eq. ( [12] ) with this expres¬ 
sion, we find the recursive relation tor the values of the 
wave function on G sites: 


fjCj — Xj+iCj+i + XjCj 
fj = A Cj-Xj~Xj+1, Xj = 


- 1 J 


A G 2 

A Bj ~ 


±Aj 


(14) 

(15) 


For the doublet states and its Hermitian counterpart 
shown in Fig. [5j A zj = A^ since it is the same for 
all Zj lattices on the left (right) side of x = 0, and 


A^ = A c = A. We can then rewrite Eq. (14) as 


a 


■7 + 1 = 

Cj 


a b a -g 2 

J 2 


-2 


Cj-l 
Cj ' 


(16) 
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Defining the localization length as £ _1 = In \Cj/Cj^% \ « 
\n\Cj-i/Cj\ for x > 0 and assuming |A#|, |A| <C G, J, 
we find 


r 1 w In (f - V+2 - l) , 




(17) 


APPENDIX: LOCALIZATION LENGTH OF THE 
DEFECT STATES 

In the main text we estimated the localization length 
of the doublet states by a degenerate perturbation theory. 
This localization length can also be estimated, for example, 
using the iterative relation for the values of the wave 
function on C sites. To derive this iterative relation, we 
apply the effective Hamiltonian H = Hq + H' to all the 
lattice sites: 

A Aj Aj=GBj, (10) 

A BjBj = GAj + J(Cj 4 Cj- 1), (11) 

A CjCj = J(Bj + 1), (12) 


using Eq. (16), and it gives £ ~ 1.04 when G = J. This 
result holds whenever there is a defect state in the Lieb 
lattice, and it applies to the point defect in both Hi of the 
“V” configuration and the Hermitian defect we discussed 
in Fig. [5] The difference of these two cases lies in the 
values of A#, A, which are not important for £ as long 
as | A# |, | A| <C G, J. Finally, we note that by defining a 
complex wave vector k = — i£ -1 for x > 0, Eq. (17) can be 
interpreted as solving the band structure of the dispersive 
bands inversely m ■ instead of finding the energy of a 
band at a given wave vector, one can find the wave vector 
at a given energy; if this energy is in the band gap, one 
finds that k has to be purely imaginary and its absolute 
value is given by Eq. (17) when |A#|, |A| <C G, J. 
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